home *** CD-ROM | disk | FTP | other *** search
/ Magnum One / Magnum One (Mid-American Digital) (Disc Manufacturing).iso / d18 / nrpas13.arc / CHEBPC.PAS < prev    next >
Pascal/Delphi Source File  |  1991-05-01  |  654b  |  29 lines
  1. PROCEDURE chebpc(c: glcarray;VAR d: glcarray; n: integer);
  2. (* Programs using routine CHEBPC must define the type
  3. glcarray as in routine CHEBFT. *)
  4. VAR
  5.    k,j: integer;
  6.    sv: real;
  7.    dd: glcarray;
  8. BEGIN
  9.    FOR j := 1 TO n DO BEGIN
  10.       d[j] := 0.0;
  11.       dd[j] := 0.0
  12.    END;
  13.    d[1] := c[n];
  14.    FOR j := n-1 DOWNTO 2 DO BEGIN
  15.       FOR k := n-j+1 DOWNTO 2 DO BEGIN
  16.          sv := d[k];
  17.          d[k] := 2.0*d[k-1]-dd[k];
  18.          dd[k] := sv
  19.       END;
  20.       sv := d[1];
  21.       d[1] := -dd[1]+c[j];
  22.       dd[1] := sv
  23.    END;
  24.    FOR j := n DOWNTO 2 DO BEGIN
  25.       d[j] := d[j-1]-dd[j]
  26.    END;
  27.    d[1] := -dd[1]+0.5*c[1]
  28. END;
  29.